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We examine the Meissner state nonlinear electrodynamic effects on the field and angular depen- 
^^ . dence of the low temperature penetration depth, A, of superconductors in several kinds of unconven- 

tional pairing states, with nodes or deep minima ( "quasinodes" ) in the energy gap. Our calculations 
r^ , are prompted by the fact that, for typical unconventional superconducting material parameters, the 

Q ■ predicted size of these effects for A exceeds the available experimental precision for this quantity by 

O ' a much larger factor than for others. We obtain expressions for the nonlinear component of the pen- 

; I , etration depth, AA, for different two- and three- dimensional nodal or quasinodal structures. Each 

^ l' case has a characteristic signature as to its dependence on the size and orientation of the applied 

^ ' magnetic field. This shows that AA measurements can be used to elucidate the nodal or quasinodal 

^ . structure of the energy gap. For nodal lines we find that AA is linear in the applied field, while the 

dependence is quadratic for point nodes. For layered materials with YBa2Cu307_i (YBCO) type 
anisotropy, our results for the angular dependence of AA differ greatly from those for tetragonal 
materials and are in agreement with experiment. For the two- and three- dimensional quasinodal 
cases, AA is no longer proportional to a power of the field and the field and angular dependences 
^ , are not separable, with a suppression of the overall signal as the node is filled in. 

Q ■ 72.40.Hi,74.25.Nf,74.20.De 



I. INTRODUCTION 



lO . Until about fifteen years ago, the question of the determination of pairing states in superconductors was one of 

P^ ' small and purely theoretical interest, since no existing superconductors were commonly suspecteda to be in a state 

^^ , other than the standard s-wave. Only liquid '^He was knowro to exhibit p-wave pairing in its several superfluid phases. 

Since then, the situation has dramatically chapged. First, extensive studies in high temperature oxide superconductors 

(HTSC's) have led to the widespread belieia that in most cases the order parameter for these materials is at least 

"jrt ' predominantly d-wave, with lines of nodes in a quasi two-dimensional Fermi surface (FS). Whether these are true 

Cj ' nodes or very deep Muiima is, however, not really established even in the best studied compoundsu, and the situation 

I I is less clear in someO^El other cases. More recently, unconventional pairing states have been proposed, on evidence of 

'■^ ■ varying strength, for a plethora of other materials wii|li-lower suposcotiducting transition temperatures. Among these 

Ch , materials are some heavy fermion (HF) compounda3~ll£l, mcmberstinlj of certain superconducting families of organic 

O ' saitersiich as k - (BEDT - TTF)2Cu(NCS)2 and (TMTSF)2X (X = PFg, CIO4, etc.), and^yttain other salts such 

.. aalZl^Eil Sr2Ru04 for_which a pairing state similar to that in the A-phase of '^He has beencdiSj suggested, although 

K*" • recently other statesE^ have also gained favor. 

To probe the bulk order pacameter (OP) it is best to use experimental techniques that measure properties over a 
scale of the penetration depthEO A, which in most materials of interest, as mentioned above, is much larger than the 
C^ ' coherence length ^. The OP at the surface may differ from that in the bulk, and furthermore, surface experiments 
are subject to uncertainties arising from surface quality and preparation problems. In cases where the OP leads to aa. 
energy gap which has nodes (or very deep minima which we denote as "quasinodes"), it was pointed out eight yeMS — 
ago, that excitation of quasiparticles near the nodes by an applied-magnetic field leads to nonlinear anisotropiesE^L^ 
in the electromagnetic properties of the material. It was shownEfl that one can in principle use these anisotropics 
to perform "node spectroscopy", that is, not only to detect the existence of nodes (or quasinodes) but to infer their 
location on the FS. r--|r— , 

Although in the earliest workE3'E2l anisotropics in the penetration depth were considered possible subjects of experi- 
mental investigation, emphasis soon switched to related quantities, chiefly the anisotropic component of the magnetic 



moment transverse to the applied field, and the torque associated with it. This quantity seemed more accessible 
experimentally than A, which was deemed to be more difficult to measure with the requisite precision. However, 
recent advanoea and refinements in experimental techniques force a reconsideration of this assessment. The best 
measurementfO of the transverse magnetic moment yield only a relatively weak lower bound on the magnitude of the 
gap at a quasinode, because the noise of the measurements is relatively high, with a resulting uncertainty a factor of 
only three below the maximum signal expected for a system with pure nodeS|. On the other hand, measurements of 
the penetration depth in compounds such as YBCO can now be performedEj with a precision of a small fraction of 
an Angstroni_As we show later, this is one or two orders of magnitude below the putative signal for that compound. 
FurthermoretZI, for many of the non-HTSC compounds mentioned above, the predicted nonlinear signal can be con- 
siderably larger than that expected for HTSC's. Hence, extension of techniques such as those of Ref. |3g to dilution 
refrigeration temperatures is likely to allow the determination of the nodal structure of these materials. It appears, 
therefore, that measurements of A arc the most promising way of probing the nodal structure of the gap through 
nonlinear electromagnetic effects. 

With this in mind, we discuss here the nodal spectroscopy of the penetration depth. We define the quantity 
A\(ip) = X{ip, Ha) — A('0, 0) where Ha is the magnitude of the applied field and ip the angle it forms with a suitably 
defined axis. In defining \(tP,Ha) attention must be paid to the experimental methods involved. We then present 
results that show that the field and angular dependence of AA('0) will reflect the nodal or quasinodal structure of 
the energy gap. We first consider quasi two-dimensional sytems with d-wave pairing. In this case, we show that the 
orthorhombicity as it occurs in compounds of the YBCO type, must be included in a proper and clear way. We calculate 
the necessary fields and find that the effects of orthorhombicity in AA('!/') are very important, and that their neglect 
has led to misleading conclusions in interpreting experiments. We then turn to other quasi two-dimensional nodal aiid. 
quasinodal cases, for which the required results are readily obtained from previously published field distributions.cl 
Finally, we consider three dimensional systems with point or line nodalE3 structures with quasinodal admixtures. In 
our conclusions, we point out that the penetration depth can indeed be used to perform node spectroscopy in all these 
cases and elaborate on the use of our results to interpret existing or future experiments. 

II. METHODS 
A. Penetration depth 

As explained above, our focus is on the angular dependence of the penetration depth A, as a very powerful probe 
in the understanding of the symmetry of .thp bulk OP for unconventional superconductors. The presence of nodes or 
quasinodes in the energy gap gives risec3Ej to nonlinear corrections in the current response to an applied magnetic 
field. Ha. These nonlinear corrections result in A having an angular and field dependence that reflects directly the 
symmetry of the pairing state. 

The first question we must address is that of defining the angular dependent penetration depth in the nonlinear 
case. This involves both theoretical and experimental difficulties. Let us consider the geometry of a semi-infinite 
superconducting slab of thickness d, much larger than any relevant penetration depth. The slab is oriented perpen- 
dicular to one of the symmetry axes, e.g., the c axis, and assuming orthorhombic or higher symmetry, its surfaces are 
parallel to the plane spanned by the other two (e.g. the a and h) axes. This is the geometry that we will consider in 
this work. In the linear case, the penetration depth is described in termg23 of the superfluid density tensor, and its 
principal axes are those of symmetry. The principal values Aq and Ab can be determined by experimentSriavolving an 
applied field along the h and a directions, respectively. Thus one can use, for example, as a conventional definition, 
Ai = J^ dzHj{z)/H{0), where Hj{z) is the magnetic field along a principal axis, which depends exponentially on 
the distance z within the sample through only one of the principal values. Alternatively, one can write \i in terms 
of the spatial derivatives of the fields at the surface, since both the spatial extent to which Ha penetrates into the 
superconductor and its derivative at the surface are determined by the same length. Because of the tensorial nature 
of the penetration depth, once the principal values are determined, the results of any experiment involving applying 
Ha along an arbitrary angle V' with respect to e.g., the a axis can be elucidated, although the result will depend on 
the specific experiment considered. 

In the presence of nonlinear effects, the situation is much more corpaiicated. Because of the nonlinearities, the 
superfiuid density or the penetration length are no longer tensors. Also&irEZi the lengths that characterize the surface 
derivatives and the extent of field penetration differ by significant numerical factors. Several possible definitions of the 
effective A(V') which coincide in the linear limit give different results for \[i}),Ha) in the nonlinear case. In general, 
experimentally one measures the extent to which fields penetrate and definitions involving surface derivatives are not 
appropriate. To find the right definition one must consider the experimental setup. In the experiments of Ref. |3q. 



the crystal is rotated to different orientations with respect to the field and a measurement of the compO|aept of the 
magnetic moment m along the applied field is performed. The penetration depth X{ip) is then extractecoH through 
the relation: 

II HAVf 2A(^)\ 

«"W = "^;^(1-^-) (2.1) 



where m" is the component of m along the field and V the volume of the sample. The term in 2.1 involving the effective 
penetration depth depends only on the sample area A in the direction parallel to the field. One can equivalently write 
this definitinn.pf. A(^, Ha) in terms of the integral of the appropriate component of the field by making use of standard 
identities Bj'OEJ One then has 

fd/2 

\{i^,HA)= dzH\i,,HA,z)/H{{)), (2.2) 

Jo 

where i?" is the component of H(z) parallel to the applied field. This is the definition of \{tp) we will use. Other 
definitions may have to be employed for different experimental setups. However, once the nonlinear field distributions 
inside the sample are known it is a rather easy matter, as seen below, to extract the effective A corresponding to any 
other alternative definition. 

To separate the nonlinear effects we write: 

\{^,HA)^Xun{-4>) + /^X{-^,HA) (2.3) 

Since the linear part is field independent while the nonlinear part, as we shall see below, vanishes at zero field, one 
has that the nonhnear part AX{tp,HA) = X{iP,Ha) — X{iP,Ha = 0). We will see that Xun has the expected angular 
behavior. The nonlinear part can be written in terms of the nonlinear magnetic moment runi as: 

AA(^, Ha) = -^^l (^, Ha), (2.4) 

where ml , is the parallel component of nini . Since nini is an extensive quantity, proportional to the sample area, we 
see from ( ^.4[ ) that AA is intensive. Usipg the expression for the magnetic moment in terms of the current field d, 
m = ^ J dr(r x j), standard identitiesEJO and the London and Maxwell equations, (2^) can be expressed entirelyO 



in terms of the value of the nonlinear flow field v at the surface of the sample. This will be done explicitly for several 
cases. We first however, give a brief outline below on the procedure for calculating the nonlinear fields in the cases 
where they are not yet known. 

B. Calculation of the fields 

In some of the cases of interest we will be able to use the field distributions found in previous worUSil, but in several 
others the fields must be calculated. We explain here briefly the method involved, with details of the calculations in 
the Appendices. ^-. i— . 

Within the frameworkcfl of the nonlinear Meissner effect, the relationEd between j and v is a sum or linear and 
nonlinear parts, j(v) = JHn(v) +jn; (v), where v is the flow fleld, and j is the supercurrent__.Tlie linear part is the usual 
relation jun = —epv, where p is the superfluid density tensor, while the nonlinear termEj'Eil in the low temperature 
limit considered here is: 



j„;(v) = -2eN] / dh n(s)v/V'(v;.v)2-|A(,s)|2e(-v/ • v - |A(,s)|). (2.5) 

Jfs 

Here N't is the total density of states at the Fermi level, and n{s) is the local density of states at the point s on the 



Fermi surface (FS), normalized to unity. The step function in (2.5) restricts the integration over the FS by 

|A(5)|+v;.v<0. (2.6) 

The functional relationship between j and v is then combined with the Maxwell-London equation,c3£j 

VxVxv-^j(v). (2.7) 

Using the boundary conditions, V x v|(i/2 = (e/c)HA, and v(0) = 0, we can, for the geometry under consideration. 



solve (2.7) analytically for the necessary fields and then extract AA from (2.4) 



III. RESULTS 

It is convenient to introduce several dinicnsionlcss quantities used in the calculations. These are the dimensionless 
coordinate C, which represents the coordinate perpendicular to the plane surface of the sample (measured from its 
midpoint) in units of Ai, the appropriate penetration depth tensor component in the sample plane. Its surface value 
is (s,i = d/{2Xi). The dimensionless nonhnear flow field, 

Unl,i = — ^, (o.lj 

V 

is normalized by the characteristic linear velocity, 

V = -XHa, (3.2) 

c 

where A is defined in each case from the in-plane components of the linear penetration depth. 

A. 2-D nodal lines with YBCO type othorhombicity 

As our first example, we examine the nonlinear|-eSects associated with a two-dimensional gap that has nodal lines, 
with crystal orthorhombicity of the YBCO type,EjE3 that is, with the nonequivalent a and b axes being along the 
antinodal directions. The applied field is in the a — b plane, forming an angle ip with the a axis, so that the fields have 
a and b components, which depend only on t he coordinate z. Considering the usual linear term only, it is completely 



elementary to verify that the definition (2.1) yields Xiin — Af, cos^ ip + Xa sin ijj, independent of the magnitude Ha- 
We can turn then to the nonlinear AA. 

The four line nodes are symmetrically placed at angles (pn (measured from the positive a-axis), where n = 1,...,4 
labels the node. The Fermi velocity at ipi forms an angle a with the +a axis. These angles are shown in Fig. n^. In 
the presence of orthorhombic distortion, ipi need not equal 7r/4, and a does not have to be equal to ipi. One often 
characterizes the deviation of ipi from n/A by describing the order parameter as being_''d -I- s" writing for example 
A{ip) = AdC0s(2(/3) -I- As, and then introducing a separate angular variable (see Ref. ^4|), for the orientation of the 
Fermi velocity. This may be misleading, however, since it is in this case no longer accurate to classify the OP in 
terms of angular momentum waves. More important, since the nonlinear results depend-jpnly on the local properties 
at the nodal positions, there is only one relevant angular variable, which is a, regardlessc3 of which origin one wishes 
to ascribe to itJlj 

We can see from the Fig. |l| that the magnitude of v/ is the same at all nodes, and that we can restrict the angle 
Ha makes with the a axis to ip £ [0, 7r/2]. We characterize the anisotropy of the linear penetration depth tensor by 
Aq = Xa/Xb, and it suffices to take Aa > 1. Only the local properties of the OP near the nodes contribute to the 
nonlinear current, hence we express the OP near the nodes as, 

A{ip)^2Ao{ip-ipn), (3.3) 

where Aq is half of the slope of the OP near the nodes, and should not necessarily be identified with the gap maximum. 



When ip G [0,7r/2], we have (see 2.6) the possibility of quasiparticle activation (QPA) at the nodes labeled (see 
Fig. Il|) (1) and (2), or at (2) and (3). The specific nodal pair that is activated depends on which nodal v/ satisfy the 
restriction in ( |2.6| ). Anisotropy in the penetration depth tensor leads to v twisting (with increased depth from the 
surface) towards the axis with the larger penetration depth, which is a. This effect occurs at linear order and it is 
therefore very significant. The nodes which contribute to the nonlinear term, therefore, depend on the dimensionless 
coordinate Q = z/Xi within the sample. There are three cases to consider: The first is when ip S [O,-0i], where f/'i 
is the maximum tjj that will result in QPA solely at (1) and (2). This angle is very small except when A^ > 1. The 
second case is when tjj G [■0i,'02], where V'2 is the maximum ip that will give QPA at nodes (1) and (2) until a depth 
Q is reached, then there is a crossover, and subsequent QPA at nodes (2) and (3). The third region has ip € [i/'2, 7r/2], 
where the only nodes activated at any depth are at (2) and (3). 

The nonlinear curre nt-fl ow relation for the case when there is QPA at nodes (1) and (2 ) is c alculated by inserting 
the OP, Eq. ([3.3|) into (2.5). The main steps are carried out in Appendix H. We find in Eq. (A3a) for the a-component, 



Jni,a = —2ep cos asm a. (3.4) 



Here p = ^Nfv'i, is the local value of the superfluid density at the nodes, Nf is the local density of states, Vf is the 



Fermi speed at the nodes, and the local crit ical v elocity is Vc 
the a-component of the current is given by ( A4a ), 



Aq/w/. Similarly, when there is QPA at (2) and (3), 



r r 9 9 9 • 9 1 

Jni,a = — COS a [v^ cos^ " + «& sm a\ 



(3.5) 



Anal ogou s expressions can be written for the 6-component. To find the nonlinear flow field we insert the current (3.4) 
into (2.7), for Ca S [(*, (s,a] or (3.5) for (a G [0, (*]■ The solution is found pcrturbatively (to first order) using the 
previously stated boundary conditions plus continuity of the flow field, magnetic field, and current at the crossover 
point Q. The details are given in Ax)pendix ^, the main results are (B4) and-. JB10[ ). We choose A = \^Xa^b as the 
normalization in (B.2), and we takecj A = A„, its nodal value. Hq is the usualEifEj characteristic field. 



^0 = 



cAq 
eXvf' 



(3.6) 



We can now achieve our objective, and get AA from the calculated fields. We write (2.4) in terms of the dimensionless 
flow field. 



AA = A [sinipUnLaiCs.a) ~ COS 'ip Uni^b{Cs,b)] ■ 



(3.7) 



The flow fleld results in Appendix M arc valid for any material thickness d. Upon taking the slab limit d ^ A, we can 
express AA(V') in the following form: 



Hq 



(3.8) 



We see from this result that the nonlinear effect in the penetration depth is proportional to the field, for line nodes. 
The quantity y{tp) represents its angular dependence, normalized so that its maximum is unity for the tetragonal 
case. The function y has a different expression for each region of -0, with the crossover angles tpi, V'2 being given in 
(Bl). These expressions are: 



yw 
yw 



^ cos^ a sin a sin^ th cos iIj + -rrr sin"^ a cos'^ ih, th & [0, -01 

2 + Aa Kl V I , Y- 

18Aa 2 • ; ■ 2 ; 

■ COS a sm a cos ip sm jp 



2 + A„, 



(3.9a) 
(3.9b) 



sin a cos ■0 



A2(1 + 2A,) 

2A^(2A^ - lOAg - 1 
(2 + A„)(l + 2A„) 



1 + 2A, + (4A„ - 1) 



Aatan^A '^ 



tana 



3 • 3 , /'AatanV'\''" 

cos a sm -0 

tana 



■0 G [V'l,V'2], 



y{^) = 



— ;— sin^ a cos a cos^ li sin li + 2A^ cos^ as\v? ib, ^ € \tp2, ttI. 



(3.9c) 



When Aa = 1, "01 = "02 = Q^ and the middle one of the expressions above is not needed. We present plots for yiip) in 
the next three Figures, where ip is limited to < < 7r/2, since the result in the remaining range is trivially obtained 
by symmetry. In Fig. S, we show results for fixed a = 7r/4, and vary the anisotropy parameter A^ of the penetration 
depth tensor, which influences the result through the twisting of the flelds inside the sample. For reference, we plot the 
isotropic result as the b^d curve. The ratio of \/2/2 between the minima and maxima of this reference curve has been 
known for a long time.ES We consider the range 1.0 — 1.5 for A^, in increments of 0.1. We see that as this parameter 
increases, the symmetry of the curve changes, to reflect the mixing of tt/2 and n symmetries. The signal increase with 
Aq when -0 = 7r/2, can readily be understood physically, since the volume occupied by the currents will be in this 
case determined by the larger of the linear components. We can see from the Figure the signal characteristics change 
considerably with increasing anisotropy. The maxima and minima of y are no longer separated by the factor of \/2/2 
for Aa > 1.0. With sufficiently large Aa, as seen in the Figure, the signal at = 7r/4 becomes approximately the 
average of y{ip = 0) and y{ip = 7r/2)[_jrhis illustrates the high sensitivity of the penetration depth to anisotropy, as 
compared with the previously studiecEj transverse magnetic moment. The reason for this difference in the transverse 
and longitudinal behaviors is explained in the last Section. 



Next, in Figure ^, we examine the effects of varying a wliile keeping A^ — 1.0 fixed. Again, the bold curve is the 
a = 7r/4 result, and the other curves are for values a = (7r/4) ± nSa with n — 1,2,3 and Sa = 7r/80. It is seen that, 
at ?/; = 0, the signal increases as a increases above 7r/4, reflecting the increase in |v • v/|. For smaller values of a 
the effect reverses, with the curves corresponding to the same n and opposite sign being symmetric with respect to 
exchange of the a and b axes. This behavior is very sensitive to small changes in a. The effects of increasing a are in 
a sense opposite to those of increasing Aq. It is interesting therefore to see how these two effects combine. Thus, in 
Fig. |[ we plot the normalized AA for Aa — 1.3, and the same values of a as in Fig. ^. The symmetry of the curves 
for the same n is now lost, because A^ 7^ 1. The overall conclusion that one can draw from these results is that in the 
presence of even relatively moderate deviations from the tetragonal "pure d-wave" situation, the appearance of the 
AA(V') vs "0 data might reflect more a tt than a 7r/2 symmetry. This has to be kept present in analyzing experimental 



results. We will return to this point in Section IV 



B. 2-D line nodes/orthorhombicity of BSCCO type 

Next we consider the nonlinear effects associated with a two-dimensional OP with line nodes in the presence of 
orthorhombic anisotropy characteristic of Bi2Sr2CaCu208+5 (BSCCO) materials. The new orthorhombic a and b axes 
of symmetry form angles of 7r/4 with the undistorted tetragonal axes. In this case, the nodes are at angles of 7r/2 
from each other and the nodal Fermi velocities are aligned with the nodal directions, which are the principal axes of 
the system. However, the two nodes on the new a axis are not equivalent to the other two on the new b axis. The 



fields for this case were calculated in Section IIIA.2 of Ref. p4. We need only to insert these fields into (3.7) and get 
(for d » Aj), 

AX{^,HA)^l^XB{i:), (3.10) 

D Hq 



where A = vAnoAib, with X„i = (27re /c )A^/w?„j, where w/„,; are the nodal Fermi velocities, Hq is defined in ( |3.6| ) 
with vf = -JvfnaVfnb, and we take Xni — vKiXt- Here the simple angular dependence is contained in the factor B^tJj), 
and is given by 

In Fig. M, we show B{ip) for varying degrees of anisotropy, A^ = 1.0 — 1.5 in increments of 1/10. 

C. 2-D quasinodal lines 

Now we turn to the situation where there are no nodes, but rather, very deep minima in the gap function (quasin- 
odes). This can be due to a small, constant iAg or iA.j_^ component in the OP, which we assume is added to a main 
Ad 2 2 component. The OP near the quasinodes can then be written in the form 

A(v3) =^A,„i„ + 2Ao((y9-v3„), (3.12) 

where 2Ao is the slope of the OP at the minima of the gap function, and A,„iii <ti Aq is the minimum value of the 
energy gap. Here we consider only the isotropic case Aa = A;, = A, a = 7r/4, which is sufhcient to illustrate the 
changes brought about by the presence of quasinodes, rather than nodes, and where we can use previously calculated 



field distributions from Section IV of Ref. |3^. We insert these fields into (2T). The resulting nonlinear AA depends 



on two variables, besides the angle ip: the ratio h = Ha/Hq of the applied field to the characteristic field Hq, defined 
in (3.6), and the ratio 

_ Amin /Ha .„ „^ 

We find after straightforward algebra, in the d^ X limit, 

AA(V,i?A) = ^A§^Q(V;,Ac). (3.14) 

The angular and field dependences now no longer factorize and are (apart from the overall factor of h) represented 
by the function of two variables Q, which is normalized to unity at V' = 0, /i = 0: 



Q{iP,k) = {cosil^' - Kf{2K + cos'ip')e {cos-ip' - k) (3.15) 

+ (sin-0' - k)^(2k + smtp')ld (sinV"' - k) , 

written for < ip' < 7r/2 (trivially extended by symmetry to the remaining range) with tp' = ip — 7r/4. The step 
functions imply that the nonlinear effects vanish for h < S, with 6 = Amin/Ao, since then the field is not sufficiently 
strong to create quasiparticles of energy larger than the minimum gap value. The behavior of Q('0, k) is plotted in 
Fig. 0, where we show its angular dependence for several values of k, with the applied field above threshold. One caH. 
see how the filling of the node produces a fast decrease in the nonlinear effect on the penetration depth, as it doesEa 
also for the transverse moment. In Fig. 0, the variation of AA/A as a function of Ha/Hq is shown for several values 
of 6, at tjj — 7r/4. The field threshold effect is clearly seen, and can be read off directly from the curves. 

D. 3-D quasinodal points and lines 

We now examine cases where the FS is three-dimensional and the nodal or quasinodal structure of the energy 
gap involves points or lines. We will consider here the same situation for which, in the limit of pure nodes, the 
transverse magnetic moment was calculated in Ref. Bn For the sake of brevity, we will compute directly the fields in 
the quasinodal case, and consider the situation where actual nodes exist as the appropriate special limit. For both 
of the cases considered below, we will assume that the slab surfaces are parallel to the a — c plane, and that Ha is 
also in this plane, with ^ defined with respect to the c axis. Then, the nonlinear fields have only x or z components 
(depending on whether we are discussing lines or point nodes respectively, see below), which depend only on the 
coordinate y normal to the sample. Again it is elementary to show that X(ijj,Ha) = Xiinitp) + AA('i/;, iJ^i), where 
^lini'ip) = ^z sin tp + Xx cos^ Ip, independent of Ha and 

AX{i>,HA) = XzSin^pu„l^ziCs,z) - XxCosijUni^x{(s,x)- (3.16) 

We will use this result with the fields calculated below. 

1. 3-D point nodes 

We first consider an OP leading to a gap with two quasinodes at the poles along the z-axis. In this configuration, 
the nonlinear fields have only a z-component, which depends on the coordinate y. By symmetry, we can restrict our 
analysis to the quasinode at 6* = 0, where 9 is the usual polar angle. We take the form of the gap near the quasinode 
to be 

|A(0)| = (|A,„i„p + |Ap0ni/2^ 0^0^ (3.17) 

where Aj, is the slope of the OP near the node, and S = Amin/Ap ^ 1. This is the generalization of the previously 
studiedtJ OP to the quasinodal case. 

The nonlinear current response as a function of the flow field is calculated in Appendix Kl and we find, for its only 
nonzero component, 

Jnl,z^^ivl~v',f'e{v^^V,), (3.18) 

'Vp 

where Vp = Ap/vf, Vg = Anun/vf, and in three dimensions, p = ■^Nfv'l in terms of local values. This result shows 
that there are no nonlinear effects present if Vz < Vs- Since the flow field decreases with distance into the sample, 
there will be a depth in the material, which, in terms of the dimensionless variable C^ = y/Xz, we denote as C^ , where 



Vz = Vs, so that nonlinear corrections are absent for the region below Q. Inserting (3.18) into (p.7[), we get an equation 
for the flow field. This equation can be solved perturbatively to first order in the small parameter {Ha/HqY , where 
i7o = cAp/eXzVf . In Appendix]^ we find, taking into account these subtleties, the solution Uni,z{Cz) for arbitrary 



thickness d. The result is given in (B14). We then have from (3.16) and (B14), after taking the limit d ^ Xz, the 
result, valid for < ip < tt/2, 

AA(V', Ha) = ]k^V{iJj, K)e{siniP - k), (3.19) 

where the function V{ip, n) is given by, 



V{i,,n) = l^c^ln /' sin^+v/sin^^ _j^\ ^ ("^i^s ^ _ 5^2 ^^^^ j ^^^^^2 ^ „ ^2_ ^3 2o) 



Here we see the mixed angular and field dependence of the result, as in the function Q in ( 3.15 ). The step function 
parameter, involving the quantity 



Aniiii I Ha 
Ap / Ho' 



(3.21) 



reflects again the requirement that the appropriate field component exceed a minimum value. In Fig. B, we plot 
Viip, k), for various values of k with the applied field above threshold. One again can see the decrease of the effect as 
the minimum gap value increases. The inset contains the field dependence of AA/Az a.t i/j — 7r/2, for 5 = — 0.04, in 
increments of 0.01. There we can see that the field dependence above threshold is no longer parabolic, as is the case 
when A,„i,i = 0. 

2. 3-D line nodes 

Finally, we consider a three dimensional FS with a quasinodal line in the x ~ y plane. The angular dependence of 
the gap near the quasinodal line a,t 9 — 7r/2 has the form, 

\A{e)\ = (|A„,i„p + |A,(| - e)\Y/^ e « |, (3.22) 

with Amin ^ Ap. The transverse magnetic moment in the Amin = limit of this OP has been previouslyEZi studied. 
For this case, the fields only have x components, and we assume tetragonal symmetry. We calculate the nonlinear 
current in Appendix H where we find, 

jni,x = (vl - V^) Q{Vx - Vs). (3.23) 

VpVx 



Agai n, th e nonlinear current vanishes unless the flow field, u, is sufficiently large; Vx > fs- Once (3.23) is inserted 



into (2/7), we get an equation for the flow field that can be solved with a procedure identical to the point node case. 



This is done in Appendix M. We find in the thick slab limit, for the only nonzero component: 

UniAO = ^ [Eie^ + E2e-< + r(C)e^ + KOe"'^] , (3.24) 

where Hq = cAp/ eXvf, C, = y/A, A is the linear penetration depth in the a — 6 plane, and the constants Ei, E2, and 
the fu nctio ns r, and s are given in Appendix H. The nonlinear correction to the penetration depth is then obtained 
from (|3l|) and ( p4| ), as: 



3 Ho 

where C{^p, n) is normalized so that its maximum at k = is unity. 



AA = \x^C{i;, K)e(cos V^ - k), (3.25) 



£(V, n) = 3/^3 tan-i ( ^^°"' ^ ""^ ) + (cos^ i, - Ak") Vcos^V^-/^^. 



(3.26) 



The overall power law behavior at Amjn = is now linear in the field, and at finite Amjn a threshold effect is found. 
We plot these results in Fig. |9| where we display C{iIj,k), for various values of k. The resulting behavior is very 
reminiscent from that found in the previous case. The inset contains the field dependence of AA/A at f/' = 0, with 
S having the same values as in Fig. g. Because of the mixed dependence of these results on the field and angular 
variables, the curves shown are not linear, except in the case S = 0, and again the threshold values can be read off 
directly from the intercepts. 



IV. DISCUSSION AND CONCLUSIONS 

In summary, prompted by recent refinements in experimental techniques that allow very high precision measure- 
ments of AA, and by the ever increasing number of candidates for unconventional superconducting states, we have 
calculated the low temperature angular and field dependence of the nonlinear component to the penetration depth 
for several different two and three dimensional energy gaps with nodal or quasinodal structures. nm 

The expected signal for the nonlinear penetration depth effect exceeds the available experimental resolutionE3'c3 of 
one tenth of an Angstrom by a considerable factor. To see this, consider for example a tet rago nal compound with 



pure d-wave pairing and a linear penetration depth A = 1400A, in the YBCO range. From (3.8), we easily get that 
the difference SX = AA^ax — AAmin, between the minimum and maximum values of X{^p) would be, at h = 0.04, 
about 3 A, a factor of thirty better than the experimental resolution. This is an order of magnitude improvement 
when compared to the corresponding estimate, under the same assumptions, for measurements of the the transverse 
magnetic moment— where one has at bestEj a factor of three. Further, taking into account the orthorhombicity of 
YBCO by settings Aq = 1.6, which increases SX (see Fig. S), we find from the same equation and at the same h, 
SX w I4A. This is a factor of 140 above the resolution achieved in Refs. p6| , [l7[ 

For gap functions with line nodes, we have found that AA at fixed angle is proptertional to the fipkii-|Our results, 
however, are obtained in the low temperature, clean liput. Both finite temperatureEd jaad impuritieg23'E3 modify this 
linear behavior at smaller fields, where nonlocal effectstJ may have a similar infiuenccEJ Since the combined outcome 
of these effects is not at present amenable to reliable computation, it is safer to perform the experiments and to 
compare with theory at the largest possible fields, where the behavior should approach linearity. Ha can be incrcjased, 
all the way to the field of first flux penetration, ff/i, taking in this full advantage of this field being in practiceEj'E^ 
much larger than the Ginzburg-Landau estimate of Hd ■ 

Keeping this in mind, let us consider our results for two-dimensional gap functions with nodes for materials with. 
YBCO-type orthorhombicity. Some experimental results for the angular dependence of X{tp) for YBCO are available,EEl 
although only for a few selected directions. We have found that the angular dependence of X{ip) is extremely sensitive 
(see Fig. 0) to small departures from unity in the anisotropy factor A^, eventually resulting in a change in the apparent 
leading symmetry behavior of AA(V'), which then looks quite different from that found in the tetragonal case. Taking 



again A^ = 1.6, forEa YBCO, we find that 3^(V' = 0) w 0.30, y{^ = 7r/2) « 1.8, while 3^(V' = ^^/^Iff 1-2, very close 
to the average in the two main axial directions. This is precisely what it is found experimentally .L3 Because this is 
so different from what happens when Xa = Ab, it was mistakenly interpreted in the experimental work as evidence 
against, instead of for, the anisotropy found there being due to the nonlinear Meissner effect. The measured field 
dependenceEa departs considerably at small fields from linearity. This is apparentljC3c3 not due to temperature effects 
alone and we believe it is very likely attributable to impurities, since the zero field temperature dependence of A of 
the sample used departs appreciably from linearity for temperatures below about 3 K. At the largest fields, the field 
dependence extrapolates to linear, with reasonable values of Hq w 9000 gauss. Thus, these experimental results are 
consisteni— as far as their field and angular dependence, with our theory. The weak temperature dependence of these 
and othercj measurements remains, howeuer, a puzzle. It cannot be ruled out, given the complications involving the 
correct treatment of impurity averagingp^ in these materials, that the temperature has a relatively weak effect in 
the samples studied. Further experimental work in the same or other materials is needed. Preliaiinary results for 
single crystals of TI-2201 show A A having a linear magnetic field dependence that is interpreted^ as agreeing with 
theoretical expectations for the nonlinear Meissner effect. 

This strong sensitivity of AA to anisotropy (either to AapOx to a, the angle that v/, at the node, makes with the 
+a axis) would not be expected from previous calculationaS^ of the transverse magnetic moment, where the effects 
of orthorhombicity were not pronounced. The reason is that the transverse moment is constrained by symmetry to 
vanish, regardless of orthorhombicity, both at -0 = and at -0 = 7r/2, plus at one point in between. This constraint 
does not exist for a longitudinal measurement. 

We have also examined gaps with two-dimensional quasinodes, and found that the field and angular dependence are 
no longer separable. The angular and field dependence of AA is governed by a term linear in the field and by a step 
function indicating that a minimum threshold field must be applied to excite quasipa rticle s above the gap minimum. 
This is multiplied by a function of and of the parameter k, which is a ratio (see ( 3.13| )) relating the value of the 



gap minimum to the a pplied field strength. The signal decreases markedly as n increases. 



In Subsection HID, we investigated three dimensional gaps with points and line quasinodes. There again the 



nonlinear contribution to the penetration depth depends on a function of angle and of a parameter k now defined in 



(3.21), a step function, and a separate factor linear in the applied field for line nodes and quadratic for points. The 
situation is similar, as far as the field dependence, to that for the two-dimensional case. The signal decreases with 
increasing k and vanishes at threshold. For example, at k = 0.6, the nonlinear signals for both points and lines drop 
to about 25% of their maximum (k = 0) values. Even with such large admixtures, however, the signal is still likely to 



be within current experimental resolution. Let us estimate the signal for an OP with three d imen sional line nodes, 



similar to that which might occur in Sr2Rup4, or certaincl heavy fermion compounds. Using ( 3.25 ), with k = 0, we 
find, S\ — AAmax = |A/i. Using publishedEJ values for Sr2Ru04 we estimate A = 2000A and a value h = 0.3 for 
Ha = Hci ■ These values give a maximum signal of, S\ w 200A. The magnitude of the signal in this case is well above 
experimental resolution, and even with relatively large admixture leading to a substantial Amin the signal would still 
be experimentally tangible. 

We have focused here on the nonlinear effects on the angular dependence of the penetration depth, and we have 
shown the strong influence that anisotropy in the principal values of the linear penetration depth and orientation 
of the Fermi velocity has on the results. The methods presented in this paper can be readily extended to other 
nodal natterns and to include the nonlinearitics in the temporal response that arise from a_tjme-dependent magnetic 
field.E3o These phenomena are currently being investigated via microwave measurements. l3 
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APPENDIX A: CURRENTS 

1. 2-D nodal lines with orthorhombicity 

For the order parameter given by (B.3), the four line nodes are symmetrically placed (see Fig. |l|) at angles ipn, 
measured from the positive a-axis, where n = l,---,4 labels the node. The Fermi velocity at node 1 is v, — 



?;/(cosq;, sina). After making the replacement N'i Jp„(Ps n{s) -^ Nf J- dip/2Tr, Eq. (2.5) gives the contribution 



j^" , when quasiparticles at ipn are activated: 



^^(v) = -2eNf r ^ v/y(2A^^J^^T2A^, (Al) 

where t^ = iy9 — (^„, and the integration is limited by (^c = |v/ • v|/(2Ao). One finds, 

^^--^^"^[v^"^-^. (A2) 

Except when v is along a nodal Fermi velocity, in general two nodes must be considered. If the nodes at ip\ and ip^ 
are activated, we can get the total nonhnear current by adding jj^^ +j„; from (A2): 



Jnia = —2ep cos asm a, (Ada) 

Vc 

jni b = sin a [w^ cos^ a + v? sin^ a] . (A3b) 

Vc 

where we have introduced the local superfiuid density, p = (l/2)Nfv'i, and critical velocity Vc ~ Ao/vf. Likewise, if 
the nodes at ip2 and ip3 are activated, we get 

jni a = — COS a \vl cos^ ct + vl sin^ a] , (A4a) 

Vc 

jnib = 2ep sm a cos a. (A4b) 

Vc 
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2. Nonlinear current for 3-D quasinodes 

a. 3-D point quasinodes 



We examine first a gap of tlie form (3.17). By symmetry, we can restrict ourselves to the node at = since the 



contribution from 9 — ii \s identical. Thus, vj ss (0,0, w/z), and the relevant region of integration is limited by 9c, 



as determined from (vj • v)^ = |A(0c)| • In performing the integral in (2.5) we again replace N% jp„(Ps n{s) by 
Nf J„ dip9d9 / At: . This yields only a z-component to the nonlinear current, 



27r 



eNfVfA, l\^r ^^^(^2 _ 02)i/20(^^ _ ^^)^ 



27r Jo Ja 



ULz^ -7^ I df I ed9i9i-9'yeiv,-vs), (A5) 

where 9^ = [(vfV,)^ - A,^i„]/A2. We get, 



J.^^{vl-v',f'e{v.-vs), (A6) 



9 

V 
P 



where Vp = Ap/vf, Vg = Amin/i'/, and, in three dimensions, p = ^Nfvl. The step function reflects that the flow field 
Vz must be sufficiently large, Vz > fs, h^ order for nonlinear effects to be present. 

6. 3-D line quasinode 



For an energy gap as given in ( 3.22 ), where the nodal line is at 6* = n/2, Vfz = over the region of integration, 
which is then limited to \9 - tt/2\ < 9,. Here {9^ - ^jlf = [(v^ • v)^ - A^J /A^ = {(v^v^ cosr,)^ - A^^J /A^, 

where v\_ = (i^^ + w^)^'^ is the projection of v on the x — y plane, and 77 the angle between V}_ and the in-plane Vf. 
In our geometry Wy = and the only component of the nonlinear contribution to the current is along x. We have. 



Jnl.x r, 

zvr 



^ / / difd9vfCosip{9l-9^)^/'^Q{v^-Vs). (A7) 



After performing the integration over 0, this leaves an integral over i^. To find the specific limits in this integral, we 
transform the integral over i^ to one over 77. Using the relation ip = P -\-r], where /? is the (fixed) angle v± makes with 
the X axis, we find, 

eNfVf f^' 
3ni,x = ^^j dT^cos{l3 + T^)[{vjv^cosT^f~Al^^]Q{v^-Vs), (A8) 

where Lps = arccos(— Amin/w/w_L). Making use of cos/3 = Vx/v_l, we get the nonlinear contribution to the current: 

jnl,x = n-Va: («! - W^) 6(u_L " Vs). (A9) 

VpVJ_ 



APPENDIX B: PERTURBATION SOLUTION 

1. 2-D YBCO-type orthorhombicity 

Here we assume that the anisotropy factor, A^ > 1, and first examine the case when ^ S [ipi, '02], where the limiting 
angles ipi are defined in the text. We will need their expressions in the slab limit, d ^ A, which are: 

1 /tanae?-"(i-A")\ 
V-i = tan-i f j , (Bla) 

1 /tana\ ,^ , , 

1^2 - tan-i f -^ J . (Bib) 
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However, unless otherwise stated, all expressions below are for arbitrary d. Without loss of generality, we give details 
of the solution for the a-component of the nonlinear current, and simply give the results for the &-component later, 
since it follows from an identical procedure. We find after inserting (A3a) into (2.7), 






-Ua - 2euaUb COS a sin a = 0, (a € [C ' Cs,a] 



dQ 



La + s[iPg cos^ a + ul siv? a] =0, Ca e [0, C] 



(B2a) 
(B2b) 



Here e — A^ „(iJ^/iJo)cosa, Aj_„ — Xi/Xn (for i — a,b ), A„^ = ^^Nfyj. We can now solve Eqs. ( B2a ) 
and ( |B2b| ) perturbatively in the small parameter e, and write Ui — u^i + ewi.i. To zeroth order, we have 

Mo, a = Aa sin?/;sech(Cs,a)sinh(^a), and uo,6 = ^^t cosV'sech(^s_f,) sinh(^f,), where A^ = Xb/Xa- The first order 
solutions satisfy the following two equations: 



j2 

— ^ " ui.a - 2uo,aUo,b COS a sin a = 0, Ca & [CCs.a] 



da 



dQ 



Ul^a + [Uo,a COS^ a + ul f, HIV? tt] =0, Ca e [0, C*] 



(B3a) 
(B3b) 



The boundary condition on the nonlinear terms is duu/dC,i\Q^ . = 0. By requiring continuit y of the flow field, current, 
and magnetic field at the point ('*, we can obtain the first order solution u„; ^ = euig to (B3a): 



^nLa 



A? „ COS a 



Ha\ 



Cig COsh(Ca) + C2q sinh(Ca) + WaiCa) COsh(Ca) + 5a(Ca) sinh(Ca) 



(B4) 



Here the constants Cu and C2i are given by Cu = w^{C*) - Wi(0) - Wi{Ct), ^21 = -5j(Cs,i) + tanh(Cs,i)[wj(0) + 
Wi{Ci) - WiiCsj) - Wi(C*)], with Wg = -y.2Ua - HaWca, Wa = -^J.iGc,a, 9a = fJ.iGs,g, and 



WcAQ 



sinh((A, + 2)C») sinh((A, - 2)C^) 

4(A, + 2) 4(A, - 2) ' 

cosh((A, + 2)G) cosh((A, - 2)Q) cosh(A, Q) 

4(A, + 2) 4(A, -2) 2A, '' 

cosh((2A, + 1)0) cosh((2A, - 1)C0 cosh(C,) 



1 3 

UziCi) = -^ cosh(3 0) - - cosh(Ci), 



Ml 



M2 



M3 



4(2A, + 1) 

3h(3 0)- 

sin 2a sin 2'ip 
2 cosh(Cs,a) cosh(Cs,b) ■ 
Aq cos^ a sin^ tp 

cosh^(Cs,o) 
Ab sin a sin ip 

cosh^ {Cs,b) 



4(2A, - 1) 



Similarly, (B3b) gives 



Ul.g = Cza COsh(Ca) + Cig sinh(Ca) + ■Wa{C,g) COSh(Ca) + 5a(Ca) sinh(Ca), 

where C^a = -Wi(0), Cu = C2i + gi{Q) - ViiQ)- Here g^ = ^2K + MsW^s.a, and 

sinh((2A, - 1)C.) , sinh((2A, + 1)C,) sinh(C,) 



W,4Q) = 



4(2A, - 1) 



4(2A, + 1) 



y,(C.) = isinh(CO- j^sinh(3C» 



(B5a) 
(B5b) 
(B5c) 
(B5d) 
(B5e) 

(B5f) 

(B5g) 

(B6) 

(B7) 
(B8) 



The matching point Q (for (s,a ^ 1) is 
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Sa ~ ^s,a 



1 



In 



Ad tan ijj 



1 — Aq V tan a 
The 6-coniponent is found by a similar procedure. We find, 



Unl,b 



a2 • Ha 

^0 



.2 ■ Ha 
Uni,b = A^„sma— - 

Cbe [0,C], 



C'lb cosh(Cf,) + 6*26 sinh(Cf,) + Wb{Cb) cosh(Cb) + gtiCb) sinh(C6) 



Csb cosh(C6) + C4b sinh(Cb) + Wb{Cb) cosh(C6) + gbiCb) sinh(Cb) 



(B9) 

(BIO) 
(Bll) 



where, Wb = M2W^c,b + fJ-sUb, 9b = -^2^^8,6 - MsH, Wb = ^J.lGcM, and g^ = -^J-iGs^b- 

The fields calc ulated a bove are for ijj g ['0i,'02]- For ■0 G [0: V'i]j we get Uni^ by simply setting the cros sover po int 
C* = in Eqs. (^,plO|). Similarly, to find the nonlinear fields for V G [V'2,7r/2], we set Q = Cs,^ in Eqs. (^,|l|). 



2. 3-D quasi-nodes 

a. 3-D point quasinodes 

For the geometry we consider. Ha is in the a — c plane, and the nonlinear fields now have only a z-component, 
which depends on the coordinate y. The flow field decreases rapidly with distance into the sample, so that there will 
be a point in the material, Q, where Uz < n, so that nonlinear corrections are absent for dist ances belo w C*. Again, 
there is no restriction on sample thickness, unless otherwise stated. Inserting the current (A6) into (2.7) gives: 



dQ 



— Uy + e [u, — K 



,3/2 



Q{uz - k) = 0, 



where e = h ^ (Ha/Hq) , and k and Hq are defined in the text. (B12a) can be written as 



dQ 



CpUz 



dQ 



-Uz = Q, 



C.e[0,C]. 



We now solve (B13a) perturbatively to first order, and write Uni,z = ewiz. We find. 



UnlACz 



Hi 
Hi 



[Di cosh(C.) + D2 sinh(C.) + fc(C.) cosh(C.) + /(C.) sinh(C.)], 



(B12a) 

(B13a) 
(B13b) 

(B14) 



where the constants Di and D2 are found by requiring continuity of the fields at C,*, and given by Di = —k{C,*), 
D2 = — /(Cs z) + tanh(^s z)[k{Q) ~ ^(Cs z)]- The functions / and k are found by elementary methods, and are given 
by, ' ' 



fiCz) 



HCz) 



1 



Bk^ 



8 V m{i^) 



■In 



2 ( m(V') sinh(Cz) + \/ m^ {^j) sinh^ {(z) - k^ 
+ [5K2sinh(G) - 2m^{i:)smh^{Cz)] [m^ W sinh^ (Cz) - K^]^^'^]e(m{ij) smh{Cz) - nY 

v2 1 m(V') cosh(Cz) + •Wm2('!/;) sinh^(^2) — K^ 



«, ,,,-(™^(V) + «;^)Mn 
8 \m('ip) 

+ cosh(Cz) {-Am^ii}) - 5k^ + m^(V') cosh(2Cz)] 

X [mH^)smh\Cz) - K^]'^^\Q(m{^)smh{Cz) ~ ^) , 



(B15a) 



(B15b) 



where mi^ili) = sin V'/ cosh C^. 2. The matching point is found to be Q, — sinh [K/m^iji)]. In a sinfilar fashion, (BlSt) 
has the first order solution 



uiz - [/(C) - fiCs.z) + tanh(C.)(fc(C) - k{Q,z))] sinh(C.). 



(B16) 
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6. 3-D line quasinode. 



Ha is again in the a — c plane, and due to the form of the gap, the nonhnear fields now have only a a:-component, 
which depends on the coordinate y. Here Ux > n, in order for nonlinear effects to be present. These effects are 
therefore absent at distances below C,* . Inserting ([A9|) into (pT^) gives the following: 



.2n3/2, 



de 



Ix ^ {Ux — K ) 6(Ma; ~ n) — ^^ 



(B17) 



where the small para meter e is e = ft = Ha/Hq. 

We can now solve (B17) perturbatively to first order. We write Ux = uqx + euix- To avoid unnecessarily tedious 
calculations, we take the slab limit, and hence the zeroth order solution, uqx is, Uqx ~ — cos ip exp (C — C^). The first 
order solution, uix, is found from dividing the slab at C* into two portions. This allows (B17) to be written 



d Ux 
~dC 



2\3/2 



dC 



0, 

Ux = 0, 






The solution to (B18a) is found by methods similar to the point node case, and is given by 

u^x = Eie^ + E2 e-< + r(C) e^^ + s(C) e"^, 



(B18a) 
(B18b) 

(B19) 



The functions r(^) and s{C) are found by elementary methods, and the constants Ei and E2 are determined from the 
boundary conditions {duix/d()\(^^ — 0, and continuity of the flow field, current, and magnetic field at the point (^*: 



E,^{s{Cs)~s{C)~r{C))e-'^^+r{Cs) 
E2 = -siCs)e-''^' + r(G) - s{C) - r{C) 



KC) 



— 3n('0)Ktan 



n{ip) 






(B20a) 
(B20b) 

(B20c) 



e (n(V')e'^ - k) 



s(C) = - 



3k3 _i / v/ri2(7/;)e^C-K2 



71 (-0) 



(B20d) 



niw) V 



?i(V') 



where n{'4>) = cos -0 exp (—^s). Similarly, (BlSb) can be solved, with the only major difference begin the boundary 
condition Uix{0) = 0. We find. 



Ulx 



E3 (e^ 



where 



i?3-(s(Cs)-.s(C))e-2C=-r(C) + r(C). 
The matching point C* is found by equating ( B18a ) and ( B18b| ) at C*, giving C* — G + ln[K/cos'0]. 



(B21) 
(B22) 



* Electronic address: khalterQphysics.umn.edu 
^ Electronic address: otvalls@tc.unin.edu 
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[II A 



FIG. 1. Variables and labels for the OP in subsection 
(4). The Fermi velocity at node (1), v/, forms an angle a with the a-axis 
A generic flow field vector v is shown for illustrative purposes. 



The four nodal lines are labeled by the numbers (1) through 
The applied field Ha is at an angle ^p to the a-axis. 



FIG. 2. Angular dependen ce of AA for the OP in subsection III A, with YBCO-type orthorhombicity. We plot here the 
function 3^('i/') (Eqs. 3.9a -3.9c), vs. the an gle i/). This function represents the angular dependence of AA normalized by the field 
dependence and numerical prefactors in (3.8). In this Figure, the Fermi velocity direction is fixed at an angle a = 7r/4. The 
bolder line shows the result for the orthorhombicity parameter A^ = Aa/Aj — 1.0 (tetragonal limit), while A^ equals 1.1, 1.2, 
1.3, 1.4 and 1.5 in the other curves (from top to bottom at i/" = 0). 



FIG. 3. Angular dependence of AA, for the same OP as in Fig. ti. Again, the function y{'4>) is plotted. Here An — 1.0 is 
fixed, while a takes the values 7r/4 ± n(7r/80) with n = 0, 1,2,3. As in the previous Figure, the effect of orthorhombicity is 
quite strong. The bold line is the tetragonal (n — 0) case. The positive values of n are (dotted curves) above the bold curve 
at small angles, and those for negative n (dashed) are below. 

FIG. 4. Angular dependence of AA for the same OP as in Figs, ti and H. Here Aa = 1.3 is fixed, while a takes the same 
values as in the previous Figure. The solid curve corresponds to a = 7r/4 and the meaning of the dotted and dashed lines is as 
in the previous Figure. The combined effects of the two anisotropy parameters are seen. 



FIG. 5. Angular dependence of A A for a 2-D OP with line nodes and BSCCO-type orthorhombicity (subsection IIIB ). The 
applied field makes an angle ^/j with the new orthorhombic a-axis. The nodal lines are now along the principal axes, which are 
rotated by an angle of 7r/4 relati ve to the previous undistorted tetragonal axes. The curves are normalized so that they represent 
the function B{4>) as defined in (3.11). They correspond, from top to bottom at i/; = to values of Aa = 1.0, 1.1, 1.2, 1.3, 1.4, 1.5. 
As in previous Figures, the tetragonal limit is plotted with a bolder line. 



FIG. 6. Angular dependence of AA for a 2-D gap with quasinodes and t etrag onal symmetry as discussed in subsection [II C 



Results are normalized so that the quantity plotted is the function Q of (3.1!:). This function is plotted vs. t he an gle tp that 



the applied field forms with the a axis. The bolder line corresponds to the usual result with nodes, («: = see (3.13)), and the 
other curves, from top to bottom, are for values of k of 0.2, 0.4, 0.6 and 0.8 respectively. 



FIG. 7. Field dependence of AA (^) for the same 2-D gap with line quasinodes considered in the previous Figure. The curves 
show the quantity AA/A (see 3.14) at ?/) = 7r/4, as a function of h for values of S (see text) equal to 0, 0.01, 0.02, 0.03, 0.04 and 
0.05 These values can be read off from the threshold field values in the curves. 



FIG. 8. Angular and field dependence of AA(-0) for a 3 -D g ap with point quasinodes. The main plot is normalized so 
that the quantity plotted is the function Pj ip, k) defined in ( 3.19| ). This quantity is plotted as a function of angle for values of 
K = 0, 0.1, 0.2, 0.3, 0.4, 0.6, 0.8, as defined in (3.21 ). The bolder curve is the n — result. The inset displays the field dependence 
of AA at i/; = 7r/2. The quantity plotted is AA/Az vs. normalized applied field, for values of S equal to 0, 0.01, 0.02, 0.03, 0.04 
and 0.05. 



FIG. 9. Angular dependence of AX (tlj) for a 3-D line quasinode. The main plot is normalized so that the quantity plotted is 
the function £("0, k.) defined in (|3.2E| ). The inset has the field dependence of AA at ■;/) = 0. All parameter values are precisely 
as in the previous Figure. 
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